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Abstract 

We investigated the vibrational properties of graphene nanoribbons by means of first-principles 
calculations on the basis of density functional theory. We confirm that the phonon modes of 
graphene nanoribbons with armchair and zigzag type edges can be interpreted as fundamental 
oscillations and their overtones. These show a characteristic dependence on the nanoribbon width. 
Furthermore, we demonstrate that a mapping of the calculated F-point phonon frequencies of 
nanoribbons onto the phonon dispersion of graphene corresponds to an "unfolding" of nanoribbons' 
Brillouin zone onto that of graphene. We consider the influence of spin states with respect to the 
phonon spectra of zigzag nanoribbons and provide comparisons of our results with past studies. 

PACS numbers: 61.48.De, 63.22.-m, 63.20.D- 63.20.dk 
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I. INTRODUCTION 



The outstanding properties of graphene and graphene-related structures of nanosize gave 
rise to extensive theoretical and experimental research during the last two decades. Along 
with the heavily studied carbon nanotubes (CNT) another quasi ID-nanostructure aroused 
special interest: Terminated stripes of graphene, so called graphene nanoribbons (GNR). 
Recent progress in preparation of single layered graphene sheets^'^'^'^ allows the fabrication 
of GNRs through lithographic techniques'"^ and possibly the verification of theoretical pre- 
dictions regarding electronic and optical properties. In the course of such investigations 
interesting magnetic properties^'^'^'^'^'^'^, quasi-relativistic behavior of electrons and the 
possibility of bandgap engineering^'^'^i^ by varying ribbon widths were shown. These 
results make GNRs seem promising for future developments in nanotechnology and nano- 
electronics. 

The propagation of valence electrons in graphene structures is accompanied by excep- 
tionally strong electron-phonon coupling^^. The investigation of the vibrational spectum in 
these materials is thus of fundamental importance for the electron transport in electronic 
devices and of great general interest. 

In this paper, we present our studies of the F-point phonons of different armchair and 
zig-zag nanoribbons, obtained through ab-initio density functional theory calculations. We 
found that it is possible to classify the F-point phonon modes of hydrogen passivated GNRs 
into fundamental oscillations, overtones and C-H vibrational modes. Fundamental oscilla- 
tions and overtones can be mapped onto the graphene phonon dispersion by unfolding the 
GNR Brillouin zone onto that of graphene. Furthermore, we discuss the dependence of GNR 
phonon frequencies on the nanoribbon width. 

II. CALCULATIONS 

Graphene nanoribbons can, at least in some cases, be understood as cut and unrolled car- 
bon nanotubes. These geometric similarities suggest that the phonon spectra of comparable 
CNTs and GNRs may be similar. On the other hand, unlike CNTs, nanoribbons possess 
edges, which have a lower coordination number, and technically, they require special treat- 
ment. A widespread method to take care of the carbon atoms at the edges in calculations is 
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FIG. 1: (color online) Structure of (a) a A^-AGNR and (b) a A^-ZGNR. In each case one dimer is 
emphasized in light grey (red). We found a lattice constant of ao = 2.4656 A for a relaxed sheet 
of graphene. The ideal lattice constants of the nanoribbons are then Cag^r — VSao = 4.27 A and 
CzGNR = 0.0 = 2.4656 A. The corresponding ribbon widths, i.e. the distance between C atoms at 
opposing edges, are 'Wagnr = — l)ao and w^zqnr = — l)ao. The relaxation results in 

nanoribbon widths slightly below (for AGNRs) or above (for ZGNRs) the values calculated with 
these equations (see Tables [J and [III)- The deviations however decrease with increasing ribbon 
width. 

to passivate the edges with atoms or molecules, mainly hydrogen. Because the diversity in 
the possibilities to cut out GNRs of a graphene sheet is larger than to "roll" it into a seamless 
cylinder, i.e carbon nanotubes, the number of edge types is very large. This leads to the 
need for a classification of those graphene nanoribbons. Some approaches use a {p, q)^ type 
classification with two characteristic integers, similar to the common classification of CNTs, 
or a classification that is based on the number of honeycombs along the ribbon widtb^. 

For the purpose of this paper entirely sufficient is the approach of Son et al.— , where 
GNR are classified by their edge type into armchair graphene nanoribbons (AGNR) and zig- 
zag nanoribbons (ZGNR), and giving the number of dimers in the unit cell (see Fig. d]). 
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The unit cell with dimers is extended periodically along the 2;-direction, resulting in an 
infinitely long strip of graphene. 

We define the width of a graphene nanoribbon as the distance between the central points 
of the outmost dimers (refer to Fig.H]). The ideal ribbon width, i.e. the width of an unrelaxed 
nanoribbon, is dependent on and given by 

w^AGNR = ^{N - l)ao (1) 
and 

WZGNR = ~^(^ ^ 1)^0 (2) 

with the graphene lattice constant ag. In our calculations, the relaxed lattice constant is 
ao = 2.4656A. The width of A^-ZGNRs with odd is equivalent to the circumference of 
an {^Y^, ^^^Y^)-CNT. Thus, these nanoribbons could be rolled into armchair nanotubes. 
For A^-ZGNR with even A^, however, there are no corresponding nanotubes. Similarly, 
the width of A^-AGNRs with odd A^ is equivalent to the chiral vector of a (^,0)-CNT, 
whereas AGNRs with even A^ do not correspond to any nanotubes. 

We used density functional theory in the local approximation form^ to calculate A^- 
AGNRs and M-ZGNRs with the number of dimers per unit cell A^=7..15 and M=4..14 
respectively. Pseudopotentials were generated with the Troullier-Martins scheme^ for the 
following valence-state configurations: C 2s^(1.49), 2p^(1.50); H ls^(1.25), where the value 
in parenthesis indicates the pseudopotential core radii in bohr. The valence electrons were 
described by a double-^ basis set plus an additional polarizing orbital. The localization of 
the basis followed the standard split scheme and was controlled by an internal SIESTA—"^ 
parameter, the energy shift, for which a value of 50meV was used. This resulted in basis 
functions with a maximal extension of 3.31 A (C) and 3.2 A (H). As SIESTA works with pe- 
riodic boundary conditions, the lattice vectors in direction perpendicular to the nanoribbon 
axis were scaled in such a way that the space between periodic images of the nanoribbons 
was at least 20 A in order to prevent interaction between them. Real space integrations 
were performed on a grid with a fineness of 0.08 A, which can represent plane waves up to 
an energy of 270 Ry. 

A minimum of 30 A;-points equally spaced along the ID Brillouin zone was used. The 
phonon calculations were performed with the method of finite differences^^. We fully relaxed 
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TABLE I: : Lattice constants c and widths w of various relaxed zigzag-edged graphene nanoribbons. 
A = '^"'^''i^" ' ' is the relative deviation of the calculated lattice constant from the ideal lattice constant. 

*--ideal 



N 


c(A) 




w (A) 


■"^ideal (A) 


2 


2.461877 


-0.15% 


2.141 


2.135 


3 


2.461083 


-0.18% 


4.286 


4.27 


4 


2.461787 


-0.15% 


6.427 


6.406 


5 


2.462023 


-0.14% 


8.566 


8.541 


6 


2.46277 


-0.11% 


10.702 


10.676 


7 


2.465000 


-0.02% 


12.838 


12.811 


8 


2.46414 


-0.06% 


14.975 


14.947 


9 


2.464095 


-0.06% 


17.113 


17.082 


10 


2.464095 


-0.06% 


19.249 


19.217 


12 


2.464195 


-0.06% 


23.523 


23.488 


14 


2.464649 


-0.04% 


27.794 


27.758 


16 


2.464730 


-0.04% 


32.067 


32.029 



'^Cideal = 2.4656 A 

''calculated by Wjdeai = — l)ao 

the atomic positions of both AGNRs and ZGNRs until the atomic forces of each atom were 
less than 0.01 eV/A and minimized the total energy as function of lattice constant (refer to 
Tables [T] and [Tl| • We used a supercell approach with a 9 x 9 x supercell and the above 
paramaters to calculate the phonon dispersion of graphene. These calculations resulted 
in a F-point frequency for the modes of 1622 cm~^. This is slightly higher than the 



24. All calculated 



experimentally obtained graphene frequency of about 1580 cm~T^ 
frequencies are therefore scaled by a constant C=0.974 to achieve better comparability with 
experimental results. 



III. RESULTS AND DISCUSSION 



Graphene nanoribbons have a large length to width ratio, which results in a quasi-lD 
crystal-like behavior and is expected to lead to confinement effects for the tt orbital electrons 
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TABLE II: : Lattice constants c and widths w of various relaxed armchair-edged graphene 
nanoribbons.A = 



^ideal 


AT 
iV 


c (A) 


A G 




«^ideal (Aj 


A 

4 


4.318 


l.lTo 


3.653 


O AO 

3.698 


5 


A O T"! r\ 

4.309 


0.9% 


/I o o r 

4.885 


4.931 


6 


4.308 


0.9% 


6.107 


6.164 


7 


4.296 


0.6% 


7.353 


7.3968 


8 


/I o r\ o 

4.293 


n r 0/' 

0.5% 


8.587 


O A A 

8.629 


9 


A Cif\ A 

4.294 


0.5% 


A O A A 

9.809 


A O A 

9.862 


10 


4.289 


0.4% 


11.056 


11.095 


11 


A O O O 

4.288 


n ACr/ 

0.470 


12.185 


12.328 


12 


4.288 


0.4% 


13.511 


13.561 


lo 


A O O T 

4.287 


0.4% 


14.753 


1 A '7A /I 

14.794 


14 


A O O /I 

4.284 


0.3%) 


1 r A o o 

15.988 


16.026 




4 984 


3% 




1 7 95Q 


16 


4.284 


0.3% 


18.454 


18.492 


17 


4.284 


3% 


19.689 


19.725 


20 


4.28 


0.2% 


23.381 


23.423 


21 


4.28 


0.2% 


24.609 


24.656 


22 


4.28 


0.2% 


25.846 


25.889 


"Cideal = V3ao 

''calculated by 


4.27 A 

Wideai = 5(Af - l)ao 









perpendicular to the ribbon axis. It is therefore justified to regard the nanoribbons as 
infinitely long in our phonon calculations. Thus, the phonon wave vector in direction of the 
ribbon axis, ky, is quasi-continuous. The ribbon edges however only allow standing waves 
perpendicular to the ribbon axis, and thus induce the boundary condition 

k_L,n ■ ^ribbon = Tl ■ TT 
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on the phonon wave f{r,t) = A ■ e^'^ "^^^ leading to a quantized wave vector 

k±,n = — n (3) 

111 

*^ ribbon 

with the order of vibration n = 0..N — 1. 

We expect therefore a vibrational behavior similar to that of an elastic sheet or a 
chain of atoms with fixed or free ends, i.e., the appearance of fundamental vibrations 
and overtones. The phonon spectrum of an A^-AGNR or A^-ZGNR should comprise of 
six fundamental modes and 3-2N-6 = 6(A^-1) overtones. Therefore, in a given phonon 
spectrum, we should be able to assign overtone modes to each fundamental mode. 

A. Armchair Nanoribbons 

Our calculations yield for each AGNR a F-point phonon spectrum consisting of 3m modes, 
with m = number of atoms per unit cell. The atomic displacements of these F-point modes 
can be classified into pure longitudinal (L), transverse (T) or out-of-plane (Z) modes. Each 
F-point phonon mode can be associated with one of three types which will be discussed 
separately: (1) C-H modes resulting from the passivation with hydrogen, (2) fundamental 
modes, or (3) overtones. 

1. C-H modes 

The four hydrogen atoms in the unit cell of AGNRs give twelve vibrational modes. These 
modes show large amplitudes of the hydrogen atoms in contrast to the almost negligibly 
small displacements of the carbon atoms. They can be grouped into 6 pairs of degenerate 
modes. We find C-H modes of different polarisations at frequencies of 750 cm~^, 850-900 
cm^^, 1100-1200 cm^^ and 3100 cm^^, which are independent of the ribbon width for all 
of our studied nanoribbons. 

2. Fundamental modes 

In any nanoribbon, a group of six modes can be found that are equivalent to the six F- 
point phonon modes of graphene with respect to the phonon eigenvectors. The two in-plane 
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optical modes, in contrast to the graphene optical modes, are not found to be degenerate: the 
in-plane transverse optical mode (TO) has a higher frequency than the inplane longitudinal 
optical mode (LO) for each of our studied AGNRs. 

The frequencies of these modes are displayed in Fig. [2|, together with the frequency of 
the experimental mode in graphene. According to Son et alr^, the nanoribbons can be 
classified into families N = 3p, N = 3p + 1 and N = 3p + 2, with p a positive integer. The 
LO-TO-splitting found for ribbons of the N = 3p family is about 29 cm"^ for the the smallest 
investigated nanoribbon and about 14 cm^^ for the largest one. For the N = 3p+l family, 
we found a splitting of 12-14 cm^^ for all investigated nanoribbons. The N = 3p + 2 family 
displays a larger splitting. It is about 46 cm^^ for the 8-AGNR and decreases with increasing 
ribbon width to a value of about 27 cm~^ for a 20-AGNR. AH these LO-TO-splittings should 
be experimentally measureable. For the LO-modes with N = 3p and N = 3p + 1 and the 
TO-modes an increase of the frequency compared to graphene is found. As can be seen, 
the LO frequencies of the {3p + 2)-nanoribbons are softened. This can be attributed to the 
small band gap in the quasi-metallic {3p + 2) nanoribbons, which is smaller than 0.294 eV 
for p > 3. This is similar to the LO phonon softening in metallic carbon nanotubes^^. We 
assume that the same effect of strong electron-phonon coupling related to a Kohn anomaly 
takes place in quasi-metallic GNR— . All modes converge towards the graphene frequency 
with increasing width. 

3. Overtones 

For each fundamental oscillation, we find (A^ — 1) overtones, where the fundamental 
displacement pattern is modified by an envelope forming a standing wave with x = 1..N 
nodes. The vibrational behavior of these modes shows similarities to elastic sheets with free 
ends. The atomic displacement can be described by an enveloping cosine function 

fn — -^n COS k±^nX = v4„ COS -- X, 

where n is the order of vibration, /c_l,„, A_l,„ and An refer to the wavenumber, the wavelength 
and the amplitude of the nth order vibration. For the wavenumbers hold the following 
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FIG. 2: LO and TO fundamental mode frequencies of armchair nanoribbons. Our calculated 

phonon frequencies were scaled by a constant factor C=0.975 to achieve a better compatibility of 
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calculations and experimental values. The dashed line indicates the experimental E2g frequency of 
pranhene. 



relations. 

k±,o = 



27r 

k±,n = ^ (4) 
nn 

= (AT - I) ao 

The nodes do not have to coincide with carbon atom positions in the unit celL Figure [3] 
shows the displacement patterns of a 7-AGNR. We characterize the phonon modes by their 
direction of vibration (transverse, longitudinal, out-of-plane) and their nature (acoustic, 
optical) as n-L/T/ZA and n-L/T/ZO, with n = number of nodes. 



B. AGNR modes in relation to graphene 

Figure [4] shows the Brillouin zones of graphene, armchair nanoribbons and zig-zag 
nanoribbons. The hexagonal structure with high symmetry points K and M represents the 
Brillouin zone of graphene with the distances TK=4:n/3ao and TM=2'7i/ \^aQ. As already 
discussed, the phonon vectors in nanoribbons are restricted by an edge-induced boundary 
condition, resulting in quantized wave numbers k±^n along the ribbon width (Eq. [5|). The 
component in axial direction however is unrestricted and not quantized. We find that the 
Brillouin zone of graphene nanoribbons consists of equally spaced discrete lines, similar 
to the Brillouin zone of carbon nanotubes^. The line spacing for armchair nanoribbons is, 
from Eq. [5|, AA;^ .„ = (jv-i)ao • translation vector of an armchair nanoribbon is given by 
Sagnk = fli + 02, i-G. the axial direction of armchair nanoribbons corresponds to the FM- 
direction in graphene. The direction perpendicular to the ribbon axis corresponds to the 
FKM-direction. The AGNR F-point overtone vibrations therefore correspond to vibrations 
in TKM direction. It should be possible to "unfold" the Brillouin zone of a nanoribbon 
onto that of a graphene sheet, where the F-point frequencies of fundamental and overtone 
modes reproduce discrete graphene modes along the TKM direction. For the overtone of 

the highest order, i.e. n = — 1, we find from Eq. [5] 

2{N-l)n 

^"■"-^ - 

271 

ao 
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FIG. 3: (a) longitudinal-optical (LO) and (b) transversal-optical (TO) fundamental and overtone 
modes at the T point of a 7-AGNR. The arrows display the displacements of the atoms in the unit 
cell. The displacement strength is normalized li> emphasize the node positions. For the n-L/TO, 
the eigenvectors of the atoms reverse re-times 0-L/TO across the ribbon width. This is further 




FIG. 4: (color online) Brillouin zones of graphene, armchair (10-AGNR) and zig-zag (10-ZGNR) 
nanoribbons. Si and 02 are the lattice vectors of graphene, 61 and &2 are the reciprocal lattice 
vectors. Note that the Brillouin zone of ZGNRs ist idelaized; in actual nanoribbons it reaches the 
K-M-K' line only in the limit of large nanoribbons, see discussion in sect. Ill C. 

As can be seen in Fig. [U IFA'MI = 27r/ao, thus it should be possible to reproduce the whole 
r/CM-dispersion of graphene by nanoribbon F-point overtones. Figure [5]shov\rs a mapping of 
the resulting pairs {k^.n^^n) of AGNR F-point phonon modes onto the phonon dispersion of 
graphene. Our calculated dispersion is in very good agreement v^^ith experimentally obtained 
results^^. As we defined the longitudinal and transverse direction with respect to the ribbon 
axis, the nth overtone of a longitudinal ribbon mode corresponds to a transverse phonon at 
wavevector fc^ „ in graphene and vice versa. 
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FIG. 5: (color online) (a) Mapping of TO (filled squares), LO (circles), ZO (filled diamonds), TA 
(open triangles), LA (filled triangles), ZA (pluses) fundamental and overtone frequencies of an 
15-AGNR onto our calculated phonon dispersion of graphene (solid lines). Longitudinal ribbon 
modes correspond to transverse graphene modes, (b) Eigenvectors of selected ribbon TO and LO 
frequencies in TKM direction. The corresponding modes are indicated in (a). 
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The overtones of one respective ribbon fundamental mode reproduce different phonon 
branches of graphene in TK- and ii'M-direction. For example, the ribbon TO frequencies 
reproduce the graphene LO branch in FiT-direction, but then switch to the graphene TA- 
branch beyond the K-point. The reason for this lies in the strict mode classification applied 
to the nanoribbons. We can also understand the branch switching from the Brillouin zone of 
graphene (see Fig. HI): going along TK and then along KM, the direction of the wave vector 
changes by 120°, if one continues to stay in the Brillouin zone of graphene. Moreover, 
it is well-known that close to the K-point the phonon modes loose their purely longitudinal 
or transverse character. Figure [5] (b) shows the phonon eigenvectors of different overtones 
of the nanoribbon modes. It is clearly seen that the overtones near the graphene i^-point 
have a mixed character with displacements in different directions [see panels 3 and 4 in 
Fig. M (b)]. In total, only small deviations between the calculated graphene dispersion 
and the zonefolded nanoribbon frequencies are found. It is expected that ribbon frequencies 
converge towards the graphene dispersion with increasing ribbon width, a result which we 
find confirmed. The root mean square deviation between graphene LO branch and their 
corresponding ribbon modes decreases from 59 cm^^ (7-AGNR) to 31.5 cm^^ (14-AGNR). 

In general, the overtones reproduce the graphene phonon dispersion fairly well. However, 
they cannot mimic the special circumstances near the i^'-point. There, the TO-mode of 
graphene displays a noticeable drop in frequency due the a Kohn anomaly^"^"^, i.e. a strong 
electron-phonon coupling. As visible in Fig. [6] (a), the nanoribbon overtones show a poorer 
reproduction of this drop of frequency in the surrounding of the graphene i^'-point. This 
is understandable because of the semiconducting nature of the investigated nanoribbons, 
which would prevent the formation of Kohn anomalies. In this case, armchair nanoribbons 
with smaller band gaps should reproduce the graphene dispersion near the i^-point better 
than AGNRs with larger gaps. Indeed, our calculations suggest that the quasi-metallic 
AGNRs of the N = 3p + 2-family, which have very small band gaps, give slightly better 
results than the ribbons of the other families. As mentioned above, the nanoribbon modes 
change their character from longitudinal to transverse direction and vice versa when crossing 
the i^-point. In particular, the displacement vectors of the ribbon LO between K and M 
show strong similarities with the LA branch in graphene. Correspondingly, the frequencies 
agree well with the LA branch (Fig. [6]a). Similarly, the nanoribbon TO modes switch in 
characteristics to TA eigenvectors when crossing the -ft'-point in M-direction (Fig. [6]b). We 
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find thus that the overtones of the nanoribbons can be well mapped onto the graphene 
dispersion when including the character of their eigenvectors. In graphene, the overtones 
correspond to different branches in the T-K and K-M parts of the Brillouin zone. Overall 
we encountered difficulties in characterizing modes close to the graphene Ji -point. The 
characteristic overbending of the graphene LO-mode due to a Kohn anomaly at F-point is 
found by zone-folding of the ribbon modes, too, for nanoribbons of sufficient width (A^ > 8). 
However, the observed overbending is considerably smaller than the one in the calculated 
graphene dispersion. 

C. Zigzag nanoribbons 

Recent studies show that the ground state of zigzag nanoribbons displays antiferromag- 
netically ordered spin states^'^'^. Calculations using spin-polarization predict the opening 
of a band gap for the otherwise metallic ZGNRs and half-metallic behavior when an electric 
field is applied due to the opposite behavior of different spin directions in electric field^^. 
We analyzed the phonon spectra of ZGNRs for effects due to the bandgap opening by gen- 
erating a pseudo potential including spin polarization with using exactly the same cutoff 
radii as the pseudo potentials we used for calculations neglecting spin-polarization. We find 
excellent agreement with the results of Son et alr^ regarding the band gap. While the 
spin-polarization effects are vital for the electronic properties, we observe only small effects 
on the vibrational frequencies. Our calculated F-point frequencies differ by just up to 8 
cm~^ from calculations neglecting spin polarization. For this reason, we did not include 
spin-polarization in the following phonon calculations. 

In order to study the vibrational behavior of zigzag nanoribbons, we carried out calcu- 
lations on A^-ZGNRs with = 2.. 14 and performed the same rescaling of the calculated 
frequencies as was done for AGNRs. A distinction of the F-point phonons in fundamen- 
tal modes, overtones and C-H-modes is performed as for the armchair nanoribbons. The 
direction perpendicular to the ribbon axis reproduced by the mapping corresponds to the 
FM-direction of graphene, see Fig.Hl The wavelengths of the vibrations of equivalent carbon 
atoms over the nanoribbon width can be described by 

An = -W^ZGNR (6) 
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By Eqs. [2|, [4] and [6l we calculate a line spacing of 

2ti 27r 



V3(iV - l)ao 



The overtone of highest order, /c_l,n-i, is then 

A;±,N-i = Afc^(iV-i: 
27r 



V3(iV- l)ao 

'^-L,N-i — /- 

V oflo 



(iV-1) 



This is equal to the graphene F-M distance, as \TM\ = Therefore, we can, in theory, 

reproduce the whole TM of the graphene dispersion by unfolding the F-point phonons of 
ZGNRs of finite width. On the other hand, we determined the wavelengths of the overtones of 
our investigated ribbons by fitting a cosine function to the respective displacement patterns 
and compared them to the theoretically expected values. For small nanoribbons, we found 
considerable deviations between the wavelength of the lattice vibration of graphene at the 
M-point and the smallest wavelength that the atomic displacements of the nanoribbon can 
describe, i.e. the wavelengths of the highest order overtones. Thus, the mapping of the 
Brillouin zone of small nanoribbons cannot reproduce the whole FM-direction, as /i;_L,N-i < 
point = However, these smallest wavelengths quickly converge towards the graphene 

M-point wavelength with increasing ribbon width. We performed mappings of the phonon 
modes of ZGNRs of various widths onto the graphene phonon dispersion in FM-direction 
(Fig. [7]). Again, unfolding the ribbon overtones onto the Brillouin zone of graphene shows 
good agreement, which improves for increasing ribbon width. The overtones of highest 
order of the optical modes of small nanoribbons appear to be clinched due to the mentioned 
deviations of fitted and theoretically predicted wavelength. The acoustical modes, however, 
display a great agreement of nanoribbon overtones and graphene dispersion. The ZGNR 
fundamental mode frequencies correspond to the six F-point frequencies of graphene. For 
ZGNRs we see, in contrast to AGNRs, a clear separation between the frequencies of in-plane 
acoustic and optical phonon modes. The in-plane acoustic modes are found in a frequency 

16 



interval of 0-1300 cm~^, the (in-plane) optical modes lie between 1300-1600 cm~^ and are 
straightforward to classify as all of them are of pure longitudinal or transverse nature. The 
displacement pattern of ZGNRs shows no apparent mixing of modes, as found for AGNRs. 
As mentioned in Sect. Ill B, the mixing of modes in AGNRs occurs due to the symmetries 
at the graphene i^-point. However, there's no comparable point in the FM-direction, which 
is reproduced by the mapping of the ZGNR Brillouin zones onto the one of graphene. 

Figure [8] shows a comparison of in-plane phonon mode frequencies of A^-ZGNRs with 
A^ = 2.. 14. As was found for AGNRs, the 0-LO of the ribbon converges towards the graphene 
LO for increasing ribbon width. A similar behavior is found for the longitudinal optical 
overtones. For the frequency of the 0-TO, a non monotonic dependence on the ribbon width 
is observed. The calculated frequencies of the transverse optical overtones of low order are 
higher than those of the 0-TO. Similarly as for the AGNRs we thus find for ZGNR a (small) 
overbending for the graphene LO mode with our zone-folding method, at least for sufficiently 
large ribbon widths. The acoustic overtones of both longitudinal and transverse nature 
display an inversely proportional width dependence. In case of transverse overtones, this 
width dependence is well described by Uac fx A^~^ Longitudinal acoustic overtones, however, 
show a considerably weakening width dependence with increasing vibrational order, as can 
be seen in Fig. [8] (b). 

Finally, we calculated the phonon dispersions over the whole Brillouin zone of various 
small nanoribbons by means of a supercell approach. We used a supercell of 9 unit cells along 
the nanoribbon axis. Figure [9] (a) shows the dispersions of A^-ZGNR with A^ = 3 — 6. As 
can be seen, the dispersions feature the characteristic fourth acoustic mode of ID structures, 
which is a rotational mode around the z-axis. The displacement pattern of this mode at the 
F-point corresponds to the displacement pattern of the mode we classify as 1-ZA. However, 
the 1-ZA in our calculations has a frequency Ci;izA=5 — 20 cm"^, which we believe results from 
the presence of the hydrogen passivation and possibly also from numerical errors. The two 
out-of-plane acoustic modes converge swiftly for increasing ribbon widths, being noticeably 
separated for the 3-ZGNR, but almost degenerate for the 6-ZGNR, closely resembling the 
ZA-mode of graphene. An interesting fact was found for the phonon modes at the X-point: 
In armchair nanotubes the phonon modes are pairwise degenerate, i.e 6(n-l) phonon modes 
in (ra, ?7,)-CNTs with odd n and all modes in (n, n)-CNTs with even n, at the X-point due to 
symmetry^. Therefore, we might exiiect 6(A^ — 1) pairwise degenerate and 6 non-degenerate 
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modes in A^-ZGNR with odd A^. Similarly, all phonon modes of A^-ZGNRs with even 
should be pairwise degenerate. However, we do not find a similar degeneracy for the zigzag 
nanoribbon dispersions we studied so far. In fact, the calculated phonon spectra of A^-ZGNR 
with odd consist solely of modes that are pairwise degenerate at the X-point. In contrast, 
the phonon modes of A^-ZGNR with even A^ are largely non-degenerate at the edge of the 
Brillouin zone. 

We want to compare our results with previous calculations by Yamada et al.—. They 
applied the MO/8-method^, which uses force fields based on the Hiickel molecular orbital 
theory, i.e. a semi-empirical approach. This approach was shown to be efficient for calcu- 
lating the vibrational properties of polycyclic aromatic hydrocarbons like graphene. In this 
approach, the topology of the structures of interest is fixed consisting of hexagons with C-C 
bond lengths of 1.39 A and C-H bond lengths of 1.048 A. As can be seen in Fig. [8] and [9] 
(b), we find general agreement, but some deviations in particular for the longitudinal modes. 
We suggest the geometry relaxation that we performed in our calculations to be responsible 
for these deviations. 



IV. CONCLUSION 



We investigated the vibrational properties of graphene nanoribbons with density func- 
tional theory. We showed that the F-point phonons of graphene nanoribbons with armchair 
and zigzag type edges can be interpreted as six fundamental oscillations and their over- 
tones, which show a characteristic nanoribbon width dependence. We demonstrated that 
the F-point phonon frequencies of nanoribbons can be mapped onto the phonon dispersion 
of graphene, i.e. to an "unfolding" of the nanoribbons' Brillouin zone onto that of graphene. 
The edge magnetization and the resulting opening of a band gap in zig-zag nanoribbons has 
only a small infiuence on the phonon spectra. The behavior of overtones and fundamental 
modes for nanoribbons of increasing width was studied and a comparison of our results for 
ZGNRs with past studies performed. 
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FIG. 6: (color online) (a) n-LO and (b) n-TO overtones of iV-AGNRs with = 5 - 15 and n = 1 
(filled circles), n = 3 (empty squares), n = 5 (filled diamonds), n = 8 (empty triangles) and n = 11 
(crosses). Solid black lines are the calculated graphene modes as in Fig. [5l 
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FIG. 7: (color online) Mapping of ribbon LO (filled squares), TO (circles), ZO (filled diamonds), 
TA (down triangles), LA (filled up triangles), ZA (crosses) fundamental and overtone frequencies 
of an 8-ZGNR onto the calculated phonon dispersion of graphene (solid lines). k±^N-i does not 
reach the graphene M-point, as for ZGNRs of small width, the wavelengths obtained by fitting the 
displacement pattern of the calculated phonons are larger than the theoretically predicted ideal 
values. 
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FIG. 8: (color online) Calculated longitudinal (a,b) and transverse (c,d) F-point frequencies of the 
ZGNR in dependence of nanoribbon width. Filled (red) squares are fundamental oscillations, empty 
symbols are overtones. Solid black lines connect overtones of equal order n for different ribbons. 
Filled grey symbols connected with dashed lines show the results of calculations of Yamada et alM 
for comparison. The experimentally determined E'sg mode frequency of 1580 cm~^ is indicated by 
thick (brown) dashed lines. 
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FIG. 9: Phonon dispersions of hydrogenated 7V-ZGNRs with N=2-6 obtained by (a) DFT and (b) 
MO/8 calculations done by Yamada et alM. The dashed lines in (a) indicate the fourth acoustic 
mode, typical for ID-crystals. Dashed lines in (b) are out-of-plane vibrations. 
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